/*******************************************************************************
																				
	DESCRIPTION:  	This file produces panel C of Figure 2.
					
*******************************************************************************/

clear all
global id_code 112

local model Full
local year 2006 
local vars _
		 	
* Load the data:
use "${data}/003_MainWithEnsemblePred_`model'`vars'`year'.dta", clear	
drop if missing(emplAft6M_0M_In) | missing(p_emplAft6M_0M_In)	

 * Produce kernel density plot:
twoway ///
	(kdensity p_emplAft3M_0M_In, fcolor(orange_red%30) lcolor(orange_red) recast(area)) ///
	(kdensity p_emplAft12M_0M_In, fcolor("61 11 55%30") lcolor("61 11 55") recast(area)) ///
	(kdensity p_emplAft6M_0M_In, recast(area)  fcolor(ebblue%30) lcolor(ebblue)), ///
	graphregion(color(white)) plotregion(margin(b=0 l=0))						///
	xtitle("Job-Finding Probability at Start of Spell (2006 Sample, Baseline Model)") ytitle("Density") 			///
	xscale(titlegap(2)) yscale(titlegap(2))	xlabel(0(0.1)1) ylabel(0(1)8, angle(0) format(%9.0g)) ///
	legend(rows(1) order(1 "3-Month" 3 "6-Month" 2 "12-Month") region(lwidth(none))) ///
	name(distFull06, replace)

graph export "${output}/${id_code}_PredictedJFR_Distribution_DifferentJobFindingHorizons_Full_2006_kdensity.pdf", as(pdf) replace

 

 